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' Abstract. The Astrodynamical Space Test of Relativity using Optical Devices (ASTROD) 

. mission consists of three spacecraft in separate solar orbits and carries out laser interferometric 

' ranging. ASTROD aims at testing relativistic gravity, measuring the solar system and 

detecting gravitational waves. Because of the larger arm length, the sensitivity of ASTROD to 
gravitational waves is estimated to be about 30 times better than Laser Interferometer Space 

(N . 

Antenna (LISA) in the frequency range lower than about 0.1 mHz. ASTROD I is a simple 

, version of ASTROD, employing one spacecraft in a solar orbit. It is the first step for ASTROD 

I ' and serves as a technology demonstration mission for ASTROD. In addition, several scientific 

^ ' results are expected in the ASTROD I experiment. The required acceleration noise level of 

, ASTROD I is 10"" m s"^ at the frequency of 0.1 mHz. In this paper, we focus on local gravity 

gradient noise that could be one of the largest acceleration disturbances in the ASTROD I 

l/^ I experiment. We have carried out gravitational modelling for the current test-mass design and 

. simplified configurations of ASTROD I by using an analytical method and the Monte Carlo 

method. Our analyses can be applied to figure out the optimal designs of the test mass and the 

^H, constructing materials of the spacecraft, and the configuration of compensation mass to reduce 

I ■ local gravity gradients. 
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■ 1. Introduction 

5^ I A gravitational mission, Astrodynamical Space Test of Relativity using Optical Devices 

(ASTROD) ^El; was proposed to test relativistic gravity, measure the solar-system parameters 
with high precision and detect gravitational waves from massive black holes and galactic binary 
stars. The concept of ASTROD is to put two spacecraft in separate solar orbits and carry 
out laser interferometic ranging with a spacecraft near the L1/L2 point. A simple version of 
ASTROD, ASTROD I, has been studied as the first step to ASTROD. ASTROD I employs one 
spacecraft in a solar orbit and carries out interferometric ranging and pulse ranging with ground 
stations [H]. 

The acceleration disturbance goal of ASTROD I is 10"^'^ m s~^ Hz~^/^ at frequency v 
of 0.1 mHz. Assuming a 10 ps timing accuracy and the acceleration noise of 10"^'^ m s~^ 
Hz~^/^ at frequency of about 0.1 mHz, a simulation for 400 days (350-750 days after launch) 
showed that ASTROD I could determine the relativistic parameters 7 and /?, and the solar 
quadrupole parameter J2 to levels of 10~^, 10~^ and 10~^, respectively jl]. In order to achieve 
the acceleration disturbance goal, a drag-free control system using capacitive sensors will be 
employed. 

A preliminary overview of sources and magnitude of acceleration disturbances for ASTROD I 



is given by Shiomi and Ni According to their estimates assuming simple models, local gravity 
gradients can be one of the largest contributions to acceleration disturbances in ASTROD I. 
Therefore, an elaborate gravitational modelling seems necessary. 

The sources of acceleration disturbances due to local gravity gradients can be classified into 
two categories. One is thermal deformation of the spacecraft and the payload, mainly caused 
by solar radiation. Inherent temperature fluctuations in solar radiation result in unwanted 
acceleration. Because composing materials of the spacecraft and payload vary in thermal 
expansion coefficient, there would be complicated relative positional changes inside of the 
spacecraft. Elaborate thermal modelling works are required for a complete analysis. However, 
we will not discuss this in this paper. The other is positional fluctuations of the test mass. Even 
when there is no deformation in the spacecraft and payload, positional fluctuations of the test 
mass produce unwanted acceleration. 

Xu and Ni did a preliminary work of gravitational modelling for the ASTROD mission |^. 
They calculated gravitational interaction between a test mass and a cylindrical reference mass (a 
hollow cylinder with two end disks) by using the expression derived for the shape design of STEP 
(Satellite Test of the Equivalence Principle [7|) test-masses ;S^. They evaluated the magnitude of 
gravitational acceleration caused when the test mass, located at the centre of the reference mass, 
was shifted along the axial axis of the reference mass. The expression they used was obtained for 
the analyses of cylindrical bodies with homogeneous density. We use more general expressions 
(section 131) , applicable to arbitrary shapes, and consider the gravitational acceleration between 
the ASTROD I test mass and cylindrical bodies (section and rectangular parallelepiped 
objects (section IH)) in this paper. Also, we carry out the Monte Carlo simulation (section 
to estimate the gravitational acceleration due to positional fluctuations of the ASTROD I test 
mass. 

2. The configuration of tlie ASTROD I spacecraft 

The ASTROD I spacecraft has a cylindrical shape with diameter 2.5 m and height 2 m. Its 
surface is covered with solar panels. The cylindrical axis is perpendicular to the orbit plane 
and a telescope on board is set to point toward a ground laser station. The total mass of the 
spacecraft is about 350 kg and that of payload is 100-120 kg (see jH] and JHl for more detailed 
descriptions of the configuration). The orbit distance from the Sun varies from about 0.5 AU to 
1 AU (see figure 2 of [Hj for a detailed description). 

The test mass {Mtm = 1-75 kg) is a rectangular parallelepiped (50 x 50 x 35 mm'^) made 
from Au-Pt alloy with density of 2 x 10^ kg m~'^. The test mass is located at the centre of the 
spacecraft. The six sides of the test mass are surrounded by electrodes mounted on the housing 
for capacitive sensing. The gap between each side of the test mass and the opposing electrode 
is 2 mm. 

3. Acceleration of a test mass 

The gravitational potential energy of a test mass (density distribution /Jt(x') and volume vt) in 
a gravitational field produced by a source mass (density distribution /9s(x) and volume Vs) can 
be written by: 




oo 




where 





(3) 



Inner gravitational multipole moments and outer gravitational multipole moments, qim and Qim, 
represent the mass distribution of the test mass and the source mass, respectively. G (= 6.67 x 
10^^^ N m^ kg~^) is the gravitational constant. 

The force between the test mass and the source mass in the sensitive axis (say, the x-axis) 
can be obtained by shifting the multipole moments of the test mass along the axis by dX' . This 
method was used by Speake to obtain the z-component of force for STEP test masses A 
detailed description of the expression for STEP is given in section 3.2 of |12j . 

Using the formula by D'Urso and Adelberger (equation (10) of |13j^. the leading order term 
of the shifted multipole moments is given by: 



1 /(2/ + 3)(Z±m + l)(Z±m + 2) , . 
for L = l+l and M = mil. With equation the x-component of the force is given by: 
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Prom equation (jlj, one can see that a positional fluctuation Xp of the test mass produces 
the first leading terms of qim, which are proportional to qooXp. From equation ©, one can see 
that these terms couple to Q2,m±i and produce unwanted acceleration. It should be noted that 
the magnitude of the unwanted acceleration is independent of the shapes of the test mass to the 
first order, but is dependent on the mass distribution of the spacecraft {Q2,m±i)- This indicates 
that the mass distribution of the spacecraft has to be designed carefully. 

We apply these formulae to estimate gravitational acceleration of the ASTROD I test mass. 



4. Monte Carlo integration 

The x-component of force of a test mass can also be obtained by calculating the following 
formula: 

F^ = g( dx'j rf,^3 /'^(^Op.(x)(x-xO 

JVs 'JVt |X X I 

We have carried out the integration over the volumes by the Monte Carlo method ^1]. A pair 
of random points, one is in the region of vt and the other is in the region of a box that includes 
Vs-, was generated at least 10^ times in each simulation. The integrand was estimated and added 
up every time when the random point generated in the box was inside of the region of Vg- The 
same simulation was repeated at least 70 times. The average and the standard deviation of the 
results were estimated. 



5. Cylindrical Spacecraft and the ASTROD I test mass 

We consider the gravitational acceleration between the ASTROD I test mass, located at the 
centre of the spacecraft, and the spacecraft. The test mass is a rectangular parallelepiped. 
Assuming a uniform density for the test mass, the following terms of gravitational multipole 
moments of the test mass (with the origin at the centre of mass of the test mass) are null because 
of the geometrical symmetries: /-odd terms, m-odd terms and terms with m = 2,6,10,14,.... 
Therefore, the leading terms of the gravitational multipole moments of the test mass are qqq, 
^20) Im-, '74,±4 and so on. When there is a positional fluctuation of Xp along the sensitive axis, goo 
and ^20 produce the leading terms gi,±i and g3,±i, respectively. From equation ©, qi,±i couple 
to Q20 and Q2,±2, and qz,±i couple to (^40 and Qi^±2- Q2,±2 and Qi^±2 are zero for homogeneous 



cylindrical bodies (with the origin at their centre of mass) because of its geometrical symmetry. 
Therefore, the acceleration of the test mass is given by 



where the relations, qn = —91,-1 and q^i = — g3,_i, are used. By substituting the following 
relations, obtained from equation @, into equation Q, 



where is a coupling constant. 

We assume that the outer dimensions of the spacecraft is 2.0 m long by a diameter of 2.5 
m and the thickness is 5 mm, and it has a uniform density of 2300 kg m^^; the mass of the 
spacecraft is 292 kg. For this spacecraft, Q20 = 3.31 kg m^^ and Q^q = 4.07 kg m^^. For the 
ASTROD I test mass, goo = = 0.493 kg and 520 = -1-17 x 10^^ kg m^. Therefore, from 

equation @, « 3.5 x 10~^^ s~^. The first term of equation ® is larger than the second 
term by more than three orders of magnitude. Therefore, in order to reduce the magnitude of 
this gravitational coupling, Q20 has to be minimized. 

To shorten the simulation time and decrease the uncertainty, we have carried out the Monte 
Carlo integration for a smaller cylinder (its inner radius, outer radius, outer length and thickness 
of each end disk are 45 mm, 50 mm, 100 mm and 26 mm, respectively). The density of the 
cylinder is assumed to be 3000 kg m~^. 10^ random points were generated in the region of the 
test mass and the region of a box (100 x 100 x 100 mm'^), in which the cylinder can fit. The 
simulation was repeated 70 times. From the average of the 70 results, we have obtained ax = 
— 2.85 X 10^^*^ m s~^ when the test mass was shifted by 1 mm along the x-axis. The standard 
deviation of the 70 results was 1.4 %. For this cylinder, Q20 = 2.80 x 10^ kg m'^ and Q^q = 
2.71 X 10^ kg m~^. Therefore, using equation ®, we obtain a-x = — 2.8 x 10"^*^ m s~^ for Xp 
= 1 mm. The ratio of the magnitude of the first and second terms in equation (jUJ is ~ 4.5 % in 
this configuration. 

6. A rectangular-parallelepiped box and the ASTROD-I test mass 

We consider the gravitational interaction between the test mass and a rectangular-parallelepiped 
box. We assume that the test mass is enclosed in the box. The gravitational acceleration between 
them has the similar relation with equation Q when the x and y dimensions of the box are the 
same; Q2,±2 and (54,±2 are also zero in this geometry. 

This configuration is similar to the electrode box for the capacitive sensing and the test mass. 
The gaps between the opposing sides of the test mass and the electrode box are 2 mm. Therefore, 
the inner dimensions of the electrode box are 54 x 54 x 39 mm^. The distance between the 
electrode box and the test mass is so close that we need to consider higher terms to estimate 
the gravitational acceleration between them. For a simplicity, we consider a box larger than the 
electrode box: the inner dimensions are 100 x 100 x 70 mm^. We assume that the thickness of 
each wall of the box is 5 mm and it has a uniform density of 10280 kg m~^. For this box, the 
acceleration of the test mass is given by 





(8) 



we obtain. 




(9) 



ax = - 




where g40=-2.63x 10"^ kg m'', 944 = -1.14x 10"'^ kg m"^, Q20 = 5.93 x 10^ kg m^^ Q40 = 
2.28 X 10^ kg g60=-4.57x 10* kg m"^ and Q64 = 2.47 x lO'^ kg m"'^. In equation (fTH)) . 
the first term is dominant. To reduce the acceleration disturbance, Q20 has to be minimized. 
We have carried out the Monte Carlo method for this configuration. We have generated 2 x 
10^ random points in the regions of the test mass and a box whose dimensions are identical to 
the outer dimensions of the box (110 x 110 x 80 mm'^). This simulation was repeated 70 times. 
From the average of the 70 results, we obtain ax = — 5.60 x 10~^^ m s~^ when the test mass 
was shifted by 1 mm along the x-axis. The standard deviation of the 70 results was 0.7 %. 

We have carried out the Monte Carlo simulation for the electrode box with the inner 
dimensions of 54 x 54 x 39 mm^. We assume that the thickness of each wall of the electrode 
box is 5 mm and the density is 10280 kg m~^. 2 x 10^ random points were generated in the 
region of the test mass and the region of a box whose dimensions are identical to the outer 
dimensions of the electrode box (64 x 64 x 49 mm^). This simulation was repeated 90 times 
and the average of the 90 results has resulted in a^; = — 3.93 x 10^^^ m s~^ when the test 
mass was shifted by 1 mm along the x-axis. The standard deviation of the 90 results was 1.9 % 
of the average. This result of the simulations indicates that the magnitude of the gravitational 
stiffness is ~ 4 x 10~* s~^. This is approximately same as the value of gravitational stiffness we 
currently use for the estimation of acceleration disturbances in ASTROD I and as the gravity 
gradient uncertainty used in current error estimates for LISA spurious accelerations J^. For 
this configuration, contributions from the leading terms in equation (|10() are comparable. To 
reduce the gravitational coupling, each moment has to be minimized. 

7. Discussion 

We have estimated gravitational acceleration between the ASTROD I test mass and cylindrical 
bodies, and rectangular parallelepiped bodies by using the analytical method and the Monte 
Carlo method. The results obtained by the two methods were consistent. 

From the estimations described above, one of the effective ways to reduce the gravitational 
couplings seems to reduce the magnitude of Q2o- This can be done by choosing the geometries 
of constructing materials of the spacecraft. For a hollow cylindrical body with an inner radius 
A, outer radius B, outer length L and thickness of the end disks T, Q2Q is zero when they satisfy 
the relation: ^ = 1 — 2^. For the spacecraft we have considered in section [21 when the thickness 
is 4 mm, the dominant term, proportional to qooQ20, becomes zero. For a box, Q20 is null when 
it is a cube. 

For constructing materials closer to the test mass, contributions from higher terms are 
comparable and they have to be reduced simultaneously. In practice, many parts of the mass 
distribution of the spacecraft are dictated by technical requirements. Therefore, the most 
practical way of reducing the magnitude of higher terms might be to have a test-mass shape 
design that minimizes the magnitude of qim- In this regard, one of the ideal shapes for the 
test mass is a sphere with homogeneous density; it only has goo- If Q20 of materials close 
to the spherical test mass is designed to be zero, the gravitational stiffness would be reduced 
significantly. Other more realistic favoured shapes for the test mass may be a cube or a cylinder 
with the aspect ratio of //a = \/3/2, where a and I are the radius and half-length of the cylinder; 
their 920 is zero. 

In summary, it can be said that, in order to minimize the gravitational couplings, the 
spacecraft mass distribution, Q2,m±i) and higher gravitational moments of the test mass {qim) 
have to be minimized. However, the design of the test mass and spacecraft will be a compromise 
between the requirements from the shape designs to minimize the gravitational couplings and 
other technical requirements. Therefore, an elaborate gravitational modelling for the final 
design of the test mass and spacecraft mass distribution will be necessary to ensure that the 
resultant gravitational coupling is sufficiently small. It might be impossible to model all the 



mass distributions on board precisely; complicated shaped objects are difficult to model and 
some approximation might be necessary in the modelling. This limit on modelling needs to 
be considered in the estimation of uncertainty in the gravity gradients effects. If the resultant 
gravitational coupling were found too large, it needs to be adjusted by employing compensation 
mass to reduce the local gravity gradients. The analytical method provides the knowledge of 
which gravitational moments contribute significantly to the resultant gravitational coupling. It 
would help to figure out the optimal designs for the test mass and the constructing materials, 
and the configuration of mass compensation if necessary. 

In practice, even if the test mass and the spacecraft mass distribution are designed perfectly 
to make the gravitational coupling sufficiently small, some imperfections, such as machining 
tolerances and density inhomogeneities, in the test mass and the composing materials of the 
spacecraft, could cause unwanted gravitational couplings. These practical issues also have to be 
studied carefully. 

8. Conclusions 

We have carried out gravitational modelling for the current test mass design and simplified 
configurations of ASTROD I by using the analytical method and the Monte Carlo simulation. 
In order to minimize the gravitational couplings, the spacecraft mass distribution, Q2,m±i, and 
higher gravitational moments of the test mass (qim) have to be minimized. The designs of the 
test mass and spacecraft mass distribution will be a compromise between the requirements from 
the shape designs to minimize the gravitational couplings and other technical requirements. To 
ensure that the resultant gravitational coupling is sufficiently small, an elaborate gravitational 
modelling for the final design of the test mass and spacecraft mass distribution is necessary. Our 
analyses can be applied to figure out the optimal designs of the test mass and the constructing 
materials of the spacecraft, and the configuration of compensation mass to reduce local gravity 
gradients. 
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